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The Wigner crystal of composite fermions is a strongly correlated state of complex emergent parti- 
cles, and therefore its unambiguous detection would be of significant importance. Recent observation 
of optical resonances in the vicinity of filling factor v = 1/3 has been interpreted as evidence for a 
pinned Wigner crystal of composite fermions [Zhu et ai, Phys. Rev. Lett. 105, 126803 (2010)]. We 
evaluate in a microscopic theory the shear modulus and the magnetophonon and magnet oplasmon 
dispersions of the composite fermion Wigner crystal in the vicinity of filling factors 1/3, 2/5, and 
3/7. We determine the region of stability of the crystal phase, and also relate the frequency of its 
pinning mode to that of the corresponding electron crystal near integer fillings. These results are in 
good semiquantitative agreement with experiment, and therefore support the identification of the 
optical resonance as the pinning mode of the composite fermions Wigner crystal. Our calculations 
also bring out certain puzzling features, such as a relatively small melting temperature for the com- 
posite fermion Wigner crystal, and also suggest a higher asymmetry between Wigner crystals of 
composite fermion particles and holes than that observed experimentally. 

PACS numbers: 73.43. Cd, 71.10.Pm 



I. INTRODUCTION 

The fractional quantum Hall (FQHE) was discovered^ 
when experimenters were searching for the Wigner 
crystaP (WC) . The expectation was that once the kinetic 
energy is quenched by forcing all electrons into the lowest 
Landau level (LL), the interaction energy would govern 
the physics and induce a WCP However, it turns out 
that, in a large range of filling factors of the lowest LL, 
the interaction energy favors the formation of composite 
fermions (CFs) instead, which form a liquid and pro- 
duce the FQHEhJ However, when the interparticle sepa- 
ration is large compared to the size of a localized particle, 
namely the magnetic length I, the WC should be stabi- 
lized. This condition can be achieved in two ways: (i) by 
reducing the magnetic length, which can be attained by 
going to very high fields (or low filling factors) ; and (ii) by 
increasing the interparticle separation, which can be con- 
veniently accomplished by considering systems close to 
integer fillings, where, to the zeroth order approximation, 
only the electrons or holes in the topmost partially filled 
LL are relevant. It is useful to differentiate between these 
two kinds of WCs, both because they originate from dif- 
ferent interparticle interactions, and because the density 
of the relevant particles in the latter actually depends on 
the magnetic field, vanishing at the magnetic field where 
the filling is an integer. We will therefore label these crys- 
tals "type-I" and "type-II" Wigner crystals, respectively. 
Much experimental work has been performed to investi- 
gate such crystalline statesP^D Recently, microwave res- 
onance experiment^ have extensively probed regions at 
low fillings and also near integer fillings and revealed res- 
onances in the conductivity that are interpreted as the 
pinning mode of a WC. 

Interestingly, even though the formation of composite 
fermions was responsible for inhibiting the electron WC 
in a broad range of filling factors, composite fermions 



themselves can form a WC, which is called a composite 
fermion WC (CFWC). Composite fermions^ are bound 
states of electrons and quantized vortices. They experi- 
ence an effective magnetic field B* and form Landau-like 
levels called A levels (ALs). Their filling factor v* is re- 
lated to the electron filling factor by v = v* /(2pu* ± 1), 
where 2p is the number of vortices attached to compos- 
ite fermions. As for electrons, there exist two kinds of 
CFWCs: 

Type-I CFWC: It has been shown theoretically that 
at very low fillings, the CFWC has lower energy than 
the electron WC, i.e. the actual ground state is at type- 
I CFWCp^HHI Essentially, electrons capture fewer than 
the maximum number of vortices available to them, and 
use the remaining degrees of freedom to form a WC. 
The CFWC has been shown to be an excellent approx- 
imation of the exact state in numerical diagonalization 
studiesPS These studies predict that composite fermions 
with a large number of attached vortices can occur in 
type-I CFWC, and that there should be transitions, as 
a function of the filling factor, from WC of one flavor of 
composite fermions to another. 

Type-II CFWC. The CFWC can also occur at CF fill- 
ings close to v* — n, which correspond to electron fillings 
close to v = n/(2pn ± 1). The state at v* = n contains 
an integer number of fully occupied ALs, and at nearby 
fillings the system has either a small density of compos- 
ite fermions in an otherwise empty AL or a small density 
of CF holes in an otherwise full AL. These will form a 
type-II CFWCP 

We compute in this paper the static and dynamic prop- 
erties of the type-II CFWC. The primary motivation of 
our work comes from the remarkable recent experiments 
by Zhu et ali^ where a resonance was observed in a nar- 
row region around v = 1/3 and interpreted as the pin- 
ning mode of the type-II CFWC. A comparison of the 
results of our calculation with the experimental observa- 
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tions shows a good agreement, thus supporting the for- 
mation of a CFWC. 

We study the CFWC by making an exact mapping of 
the problem of composite fermions in a partially filled 
AL to that of fermions in the lowest LL interacting via 
an effective interaction. Various quantities can then be 
computed using the familiar methods developed in Refs. 
[I MESl. ( The type-II CFWC was also studied by Scarola 
in a variational study that compared the ener- 
gies of the WC, bubble crystal, or stripe phases of com- 
posite fermions. That work, however, did not compute 
the shear modulus of the CFWC or the magnetophonon 
and magnetoplasmon dispersions, and also did not focus 
on the phase diagram of the CFWC.) Our main conclu- 
sions are as follows. We estimate the ranges of stability 
for the CFWC around the filling factors v = §,§,§• The 
theoretical range of stability around v = | is consistent 
with the interpretation of microwave resonance experi- 
mental results around | as a signature of the CFWC. 5 
We predict CFWC also around | and | , except in much 
narrower ranges of filling factor. We present results for 
the dispersions of the phonons as well as magnetoplas- 
mon of the CFWC. Using a model of Chitra et. alp3 we 
argue that the frequency of the pinning resonance of the 
CFWC is expected to be of similar magnitude as that of 
the corresponding electron WC, and also determine the 
magnetic field dependence; these are also consistent with 
experiments. 

A curious observation is that the Hartree Fock treat- 
ment is much more successful for the CFWC than for the 
electron WC. For example, Maki and ZotoJ^found, with 
similar methods as those employed below, that the elec- 
tron WC in the lowest LL is stable in the range v < 0.45 
(and 1 > v > 0.55 by particle hole symmetry in the low- 
est LL), which is grossly inconsistent with experiments 
that demonstrate WC only below approximately v < 1/7. 
On the other hand, the stability range we find for the 
CFWC is quite consistent with experiment. While this 
difference may seem surprising, it is easily understood. 
The Hartree-Fock approach fails for electrons in the low- 
est LL because they form strongly correlated states of 
composite fermions, not captured by the Hartree Fock 
approach. In contrast, the interaction between composite 
fermions is such that the CFWC state is destabilized be- 
fore composite fermions can form their own strongly cor- 
related FQHE states. For example, the 4/11 FQHE state, 
which is a fractional QHE of composite fermion a 21 " 24 ! at 
the composite fermion filling factor of v* = 1 + 1/3, is 
outside the CFWC region, which persists, according to 
our calculation, only up to v* ~ 1.10. 

The remainder of the paper is organized as follows. 
We describe the mapping from composite fermions in a 
partially filled AL to fermions in the LLL. Next, we solve 
the equations of motion for the CFWC, and calculate 
the shear modulus and the dispersion relations for the 
magnetophonon and magnetoplasmon collective modes of 
the this crystal. Finally, we conclude with a discussion of 
how our results compare with experimental observations. 



II. MODEL 

Electrons in the lowest Landau level exhibit exotic 
properties, which result from the formation of composite 
fermions, bound states of electrons and an even num- 
ber of topological vortices. At the rational filling factors 
v = 2np+i tne g roun d state contains precisely n filled 
ALs of composite fermions. Of interest in this paper are 
filling factors close to these rational values, when the CF 
filling factor is given by 

v* =n±D* (1) 

In what follows below, we assume that the magnetic field 
is sufficiently strong that the system is fully spin polar- 
ized and LL mixing is negligible. Our program proceeds 
along the following steps: 

(i) We treat the n filled ALs as inert, and formulate the 
problem in terms of either solely the composite fermions 
in the (n + l) th AL (for v* = n + v*) or the missing 
composite fermions in the n th AL (for v* — n— V*). The 
physics that we wish to investigate lies fully in the dy- 
namics of these composite fermions. These are sometimes 
called CF quasiparticles and CF quasiholes, but we will 
refer to them simply as composite fermions below, while 
remembering that their filling factor is v* , and that their 
density vanishes at v = 2v ^+i (where v* =0). We do 
not include the physics of AL mixing in our work; it will 
modify the form of the inter-CF interaction, but we do 
not believe that would significantly alter the results. 

(ii) The two-body Coulomb interaction between elec- 
trons generates a complex interaction between composite 
fermions that contains two, three and higher body terms. 
We assume that the physics is dominated by the two- 
body interaction (determined in the next section). This 
should be a good approximation at sufficiently small den- 
sities of the composite fermions; because the CF filling is 
found to be small in the entire range of stability of the 
CFWC determined below, we believe that the approxi- 
mation of the neglect of three and higher body interaction 
terms in the CF Hamiltonian does not cause significant 
corrections. There is also numerical evidence that this is 
a good approximation^, but a systematic study of the 
importance of higher order interactions between particles 
as a function of system size has not been undertaken.^ 

(iii) We define the equivalent problem of fermions at 
filling factor D* in the lowest LL interacting via an appro- 
priate effective interaction. The effective interaction used 
here has been shown to reproduce the desired interaction 
pseudopotentials with extremely good accuracy. 

(iv) We calculate the properties of the W C of fe rmions 
at v* in a standard Hartree Fock approach! 16 * 17 ! In par- 
ticular, the instability of the CFWC is signaled by the 
shear modulus becoming negative. 

(v) Disorder is neglected throughout. This makes our 
results unreliable at very small values of v* , where the 
distance between the (relevant) composite fermions is so 
large that disorder will likely dominate over interaction, 
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producing a localized or a glassy phase. However, we ex- 
pect the neglect of disorder to be a good first approxima- 
tion for weak disorder and not too large a lattice spac- 
ings; here the WC order ought to persist meaningfully 
over several lattice spacings. 
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FIG. 1: Plot of the semiclassical interaction between com- 
posite fermions near u* — 1/3, both for composite fermions 
in the second A level (labeled "quasiparticle" ) and CF quasi- 
holes in the lowest A level ("quasihole"). At long distances, 
both interactions are Coulombic, but have more complicated 
forms at short distance. 



Two-body Interaction between composite 
fermions 



The CF wave functions can be obtained from the elec- 
tron wave functions by the standard mapping: 



= Pi 
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Here a is a quantum number labeling the state, v* is the 
CF filling factor, and Plll is the LLL projection oper- 
ator. The interaction between two composite fermions 
in the (n + l) th AL is determ ined by computing the CF 
pseudopotentials j 14 * 15 * 19 * 2 ^ 2 ^ i.e. the energies of the pair 
in definite relative angular momentum states, by anal- 
ogy to the electron pseudopotentials in the n th LL. The 
pseudopotential for two electrons in relative angular mo- 
mentum L is defined a:P^ 
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where <j>2' ++ is the wave function of two electrons in the 



,th 



LL in relative angular momentum L. The CF pseu- 



dopotentials are analogously defined as 
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While for the calculation of of the electron pseudopoten- 
tials V£, it is sufficient to consider only two electrons in 
the n th LL, the evaluation of the CF pseudopotentials re- 
quires a full many-particle calculation, because the Jas- 
trow factor correlates all composite fermions with one 
another (including those in lower A levels). Nonetheless, 
using the standard methods^, the relevant integrals can 
be evaluated by the Metropolis Monte Carlo method, and 
the thermodynamic limits for the CF pseudopotentials 
can be obtaine d. We refer the reader to the literature for 
further details! 14 ' 15 ' 19 ' 2 ^ 

We next define a problem of fermions confined to the 
lowest LL that have an effective interaction V (r) that 
produces the above pseudopotentials: 



<<f>°£ ++ \V cS (r)\<5>°£ ++ >= V 1 



CF,n 
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where the superscript refers to the lowest LL. There are 
many possible choices for V eS (r) which produce the same 
effective interaction.^ We find it convenient to u se the 
effective interaction given by Lee, Scarola, and Jair i 15 | 19 l 
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The asterisk on 7-* indicates that the distance is being 
measured in units of the effective magnetic length I* = 
\J hc/e\B*\. The last term captures the long distance 
limit of the interaction between composite fermions: 
( e *) 2 A = {^+l?/r*l* = {^ ri f/r*{2n + iy/% which 
is the Coulomb interaction between particles of fractional 
charge e* = e/ (2n + 1) at a distance rW^ The number of 
terms kept in the sum, and the coefficients c^, are de- 
termined to produce a good approximation for the pseu- 
dopotentials obtained from the microscopic theory. We 
use below the q given in Ref . [15] . 

In what follows, we will use /* as the unit of length and 
e 2 /I as the unit of energy. We will omit the asterisk from 
the lengths for notational convenience. We will further 
set h = 1 and c = 1. In certain places of relevance, the 
units will be explicitly restored. 



B. Semiclassical model for CFWC 

Now, we consider the dynamics of fermions in the LLL 
interacting with an effective interaction. Since we are 
dealing with fermions in the LLL, they all have the same 
kinetic energy. Therefore, the ground state of the system 
is decided entirely by interactions between the fermions. 
Since we will be working at low filling factors, we make 
the ans atz that our system is described by a crystal 
state^T3 for which we choose the hexagonal (or trian- 
gular) lattice symmetry, which it is known to give the 
lowest energy for the 2D classical WCP The localized 
single particle states are wave packets: 
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centered at hexagonal lattice sites 



l 1 V3 

it = a n H — to, — m 
I 2 2 



Eq. ( 12 ) can be evaluated analytically using Mathemat- 



(9) 



where n,m are integers and a the lattice constant. The 
Maki-ZotosP wave function is obtained by placing a 
Gaussian wave packet at each lattice site and then anti- 
symmetrizing: 



(10) 



The problem can be recast as a semiclassical crystal 
dynamics problem by calculating the expectation value 
of the effective interaction: 



<* { R 3 }({r i })|* {flj} ({r i })> 

^ f(iJy) + V 3 (R i ,R j ,R k ) + ■■■ 



(11) 



i,j,k 



where the sums on the right are infinite sums over hexag- 
onal lattice sites, i?^ = |Hj — Rj\, and V(Rij) is given 
by: 



ica, but the resulting expression is too long to reproduce 
here. In Fig. [T] we plot the semiclassical interaction given 
by Eq. ( |12[ ) using the values of Ci for composite fermions 
around v — 1/3 given by [TS]. At short range, these in- 
teractions have a complicated form as a function of the 
distance, and can even be attractive, but at long distance 
the interaction is always repulsive. The attractive por- 
tion of the interaction between the composite fermions 
might seem to suggest that the crystal is unstable to 
pairing, but the issue is somewhat subtle because the 
gain in energy due to pairing is accompanied by a com- 
peting enhancement in the interaction energy between 
the pairs. It has been shown by an explicit calculation 
reported in Ref. [TS] that for sufficiently small values of 
v* a WC state is energetically preferred over the paired 
liquid state, although the formation of bubble crystals, 
where each lattice site supports more than one compos- 
ite fcrmion, can result from an attractive interaction in 
some filling factor ranges. We will not consider the pos- 
sibility of bubble crystals below, and assume that there 
is a single CF at each lattice site. 

III. EQUATIONS OF MOTION AND SHEAR 
MODULUS 



V(R ,_ Jrdr F cff (r) e -'- 2 / 4 (/ (^-r) - M±R t] r)) 
[ ij > ~ 4sinh(i?2./4) 

(12) 

The derivation of these expressions closely follows the 
work of Maki and Zotos^, and is given in Appendix A. 
As stated above, we neglect n > 3 body terms in Eq. (111. 



The individual two body terms are interpreted as semi- 
classical energy of particles located at the lattice sites of 
a hexagonal lattice. We have calculated the energy for 
both composite fermions around each of the filling fac- 



tors v = |, | , | using the effective real space interac- 



tions discussed above. With the particular form of the 
effective interaction used (see Eq. |7]))P^, the integrals in 



Next we calculate the lattice dynamics for a crystal 
of composite fermions in the harmonic approximation.^ 
We interpret the sum J2i<j V(Rij) as the zero-point en- 
ergy of the CFWC with composite fermions located at 
the lattice sites of a hexagonal lattice. We perturb the 
particle's equilibrium location Ri by a small, time de- 
pendent quantity Ui(t): 



Ri — > -R,; 



(13) 



We expand V(\Ri + u t — Rj — Uj\) to quadratic order 
in Ui,Uj and write the equation of motion for the 
particle: 



= eB*u 
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In writing these equations, we have assumed that the dy- 
namics of composite fermions^ corresponds to particles 
of mass to* in an effective magnetic field B* = B/(2n+l). 
The CF mass to* is unrelated to the electron mass!^ 
The factor of (l*)~ 2 = eB* results from expressing the 



lengths in units of I* , while recognizing that the argu- 
ment of V(Rij) is already in units of /*. This is consis- 
tent with the expectation that the magnetic field renor- 
malization from B to B* upon the formation of com- 
posite fermions is associated with an analogous electric 
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field renormalizationPS Assuming a solution of the form 
Ui(t) = Uke~ %ut ~ * the equations of motion in the x 
and y directions take the form: 

uj 2 u x k = [iuu* + w*$ X2 /fe)] u\ + uj* c <S> xx {k)u x k (16) 



and $ Q( 3 is given by 



<M fe ) = E & MJ^ { ^ 1 - cosCfe-H,-)). (19) 



2 y 



where 



-tww* + u3S + w'S^fc)"* (17) 



e£* _ eB 

m* (2n + l)m* ' 



(18) 



Here, a and /? denote X and y components. 



Equations (16) and (17) can be diagonalized to obtain 
the dispersion w±(fe): 



4(fe) = 



i 



[$ xa (fe) + $ ra (fe)])± A /j(o;f 



•^[^(fc) + $ ra (fc) 



; a [*2v( fc )-*~( fc )*vy(*)l (2°) 



To ascertain the regions of filling factor where the CF 
crystal is stable, we need to examine the phonon disper- 
sion in the long wavelength limit , k — > 0. In the limit of 
small k we can write $ Q/ 3 a d 16 * 17 ! 



$ Q/ 3 = 



+ (C L ~ C t ) k a kfi + 8 aP C t k' 



(21) 



where Ct is the shear modulus and Cl the coefficient 
in the longitudinal mode. Both have units of energy. 
In two dimensions, it can be shown that the dynami- 
cal matrix of a two dimensional hexagonal crystal with 
central forces characterized by two parameters, here 



conveniently chosen to be Ct and CxP^ As shown in 
Rcf. [T71 the classical value of C t for particles with 
charge e* in a magnetic field B* and filling factor v* is 
given by C to = 0.09775(^*) 1 / 2 e* 2 /'*, which, when con- 
verted into our energy units of e 2 /Z, reduces to Ct — 
0.09775(i>*) 1 / 2 /(2n+ l) 5 / 2 . The parameter C Lo is given 
by C Lq = -5C to PI Non-classical properties of the CFWC 
can be captured by measuring the variation of Ct and Cl 
with respect to their classical counterparts in the region 
where the crystal is stable. Substituting into Eq. (20) 
produces: 



c4(fc) = -(uj* 2 +u* c (-+C L + C t )k 2 ) ± d-(u>? + + C L + C t )ki) 



2 C t (— + C L )fc 4 



(22) 



where we have kept terms up to quadratic order in k. 
Expanding in powers of fc, we get: 



u J 2 + (k) = (u J *f+0(k), 



(23) 



IV. RESULTS AND DISCUSSION 

We give in this section the results of our calculation, 
and, where possible, compare them to experiments. 



A. Shear modulus and the region of stability 



w 2 (fc) = C t v*k 3 + 0(k 4 ). 



(24) 



The w+(/c) mode is called the magnetoplasmon mode and 
the w 2 (k) mode is called the magnetophonon mode. For 
negative values of Ct the magnetophonon frequency be- 
comes imaginary in the long wave length limit, indicating 
an instability of the CF crystal. 



To calculate the shear modulus, we evaluate the dy- 



namical matrix $ Q/ 3 of Eq. ( 19 1 numerically for several 
different values of k and fit the result to the small k limit 
form of Eq. ( |2~l| ). Further details are given in Appendix 
B. In Fig. [2] we show our results for the CFWC moduli 
around v = 1/3,2/5,3/7. In this figure, we plot C t and 
Cl normalized with respect to the classical shear mod- 
ulus of particles with fractional charge as a function of 
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FIG. 2: Plot of the shear modulus Ct and the coefficient in the longitudinal mode Cl for the CFWC at v* = n + i>* (left panels) 
and v* = n — v* (right panels), as a function of the partial filling factor v* — \v* — n|. The top row corresponds to u* ~ 1 
(i/ ~ 1/3), the middle row to u* ~ 2 (f « 2/5), and the bottom row to v* ~ 3 («/ « 3/7). Ct and Cl are normalized by their 
classical counterparts Ct and Cl . The vertical lines indicate the filling factor where the shear modulus becomes negative and 
the CFWC becomes unstable. 
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FIG. 3: Phase diagram for the CFWC. The shaded re- 
gions show the theoretically calculated stability regions of the 
CFWC in terms of the real filling factor v (top) and the CF 
filling factor v* (bottom). 



v* . Notice that as v* approaches zero, all of the crystal 
moduli become equal to the classical crystal moduli. We 
find that at sufficiently low filling factors (approximately 
v* < 1/12), the CF moduli are indistinguishable from 
their classical crystal counterparts. The non-classical be- 
havior of the CFWC becomes manifest only at somewhat 
larger values of 9* and, indeed, makes a significant cor- 
rection to the phase boundary marked by the position 
where the shear modulus changes its sign. 

The CFWC becomes unstable when Ct becomes neg- 
ative. We depict the stability regions around the filling 
factors 5 , § , f in Fig. pi which is one of the principal 
results of our work. The predicted range near v — 1/3 
is nicely consistent with the experiment of Zhu et o/P, 
who find a range ~ 1/3 ±.015, determined by where a 
microwave resonance is resolved. This differs by an order 
magnitude from the range 1±.15 for the electron crystal.^ 
In fact, as also noted by Zhu et al., when measured in 
terms of the CF filling, the range near 1/3 is given by 
v* = 0.88 — 1.15, which is more similar to the observed 
stability region of the v = 1 WC. Overall, our calcula- 
tions provide strong support to the identification of the 
microwave resonance experiments in terms of the CFWC. 
We also predict regions of stability for the CFWC around 
v = 2/5 and 3/7, which, however, are much narrower and 
may be harder to observe, especially considering that at 
very small v* the disorder may dominate producing a 
localized glassy phase. 

We note that we mark the region of stability as the fill- 
ing v* where the shear modulus first becomes negative as 
the v* is increased from zero. According to our calcula- 
tions the shear modulus becomes positive again at some 
larger fillings, to which we do not assign any physical sig- 
nificance - we believe that those regions are outside the 
regime where our approximations are valid. It was noted 
in previous calculations^ that composite fermions can 
possibly form bubble crystal and stripe phases at larger 
values of v* . 



B. Magnetophonon and magnetoplasmon 
dispersions 

In Fig. |1J we plot the dispersions of the magnetoplas- 
mon and magnetophonon collective modes for the CFWC 
around the irreducible element of the first Brillouin zone 
(FBZ) around v — 1/3. For small wave vectors, the dis- 
persion has the form uj oc fc 3 / 2 , which is identical to the 
small k behavior of the classical WC dispersion. The 
dispersion relations for the magnetophonons and mag- 
netoplasmons for the CFWC around | and I, shown in 
Figs. [5]and[6j have nearly the same shape and appear to 
be scaled versions of the dispersion relations shown for 
the crystals around v = 1/3. 

The magnetophonon and magnetoplasmon dispersions 
of the WC obtained above are valid for an ideal system 
with no disorder. At small wave vectors, a gap opens due 
to disorder, resulting in a pinning mode that has been 
studied extensively and is the topic of the subsequent 
section. As noted in Ref. [7J it is in principle possible 
to obtain information about the small q dispersion by 
comparing the pinning mode resonances measured from 
nearly identical samples that have different metal film 
coplanar waveguide (CPW) slot widths w (See any of 
[5]-[9] for experimental details.) Measurement of the dis- 
persion at large q is complicated by the fact that light 
ideally couples only to excitations with very small wave 
vectors. A possible method for accessing larger momenta 
through light scattering would be to use gratings or to 
exploit piezoelectric coupling to surface phonons to de- 
fine a wave length, which has proven useful in the study 
of collective mode dispersion of the liquid states of com- 
posite fermionsP^ The extremal points of the dispersion 
may also provide a signature in inelastic light scattering, 
in the presence of disorder, because of the singularity in 
the density of states at the corresponding energies.^ 



C. Frequency of the pinning mode 

A number of theoretical studies have considered the 
electromagneti c respon se of a pinned Winger solid in a 
magnetic field 120123(39] an( ^ g n( ^ the cyclotron res- 
onance is shifted by an amount related to the pinning 
frequency, and at the same time also broadened. The 
actual frequencies of the pinning mode depend on the 
details of the potential and are difficult to predict in a 
quantitative manner. We use below the model of Chitra 
et alW^, where they use a Gaussian variational method 
to study the elastic Hamiltonian of a Wigner crystal to 
deduce its properties in the presence of disorder. 

We first address the relation between the pinning mode 
frequency of the CFWC at an effective filling v* to that 
of the electron WC at the same filling; we believe that 
such a relationship might be more stable against the na- 
ture of the disorder than the actual pinning frequency. 
According to Chitra et. al.'s modeP2 for the pinning 
frequencies of the WC, different disorder length scales 
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FIG. 4: Magnetoplasmon and magnetophonon dispersion of the CFWC near v = 1/3 around the irreducible element of the first 
Brillouin zone. The cyclotron frequency is given by lu^ = eB* /m* where B* is the effective magnetic field and m* is the CF 
mass; the results are quoted in units of e 2 /el. For illustration, the dispersions are given at filling factors where the CF shear 
modulus has a maximum, i.e. where the CFWC is most robust, both below and above v* — 1. The irreducible element of the 
FBZ is shown in the upper right hand figure. 
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FIG. 5: Magnetoplasmon and magnetophonon dispersion of the CFWC near v = 2/5 around the irreducible element of the first 
Brillouin zone. The cyclotron frequency is given by lu^ = eB* /m* where B* is the effective magnetic field and m* is the CF 
mass; the results are quoted in units of e 2 /el. For illustration, the dispersions are given at filling factors where the CF shear 
modulus has a maximum, i.e. where the CFWC is most robust, both below and above v* = 2. 
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FIG. 6: Magnetoplasmon and magnetophonon dispersion of the CFWC near v = 3/7 around the irreducible element of the first 
Brillouin zone. The cyclotron frequency is given by cj^ = eB* /m* where B* is the effective magnetic field and m* is the CF 
mass; the results are quoted in units of e 2 /el. For illustration, the dispersions are given at filling factors where the CF shear 
modulus has a maximum, i.e. where the CFWC is most robust, both below and above v* = 3. 



produce significantly different B dependence for the the 
pinning frequencies. Their approximate expression for 
pinning mode frequencies is given by: 



(25) 

2 p 3 / 2 is the classi- 



where £ = c^tt 2 )" 1 / 6 i?" 2 (^) 6 , c 
cal shear modulus for the electron WC in units of energy 



density, R a 



a is the lattice spacing, A is a mea- 



sure of the electrostatic coupling of the particles in the 
lattice to disorder (and is proportional to the square of 
the charge of the particle in the lattice), £q =max[r 2 , I 2 ], 
p m = m/ita 2 is the mass density, p = (im 2 )^ 1 is the 
particle density and rj is the disorder correlation length. 
The pinning frequencies have a different B-ficld depen- 
dence depending on which of rj or / is larger. Putting 
everything together, we find 

a 4 p 2 A 



(26) 



Noting that a, p, muj c , and £o ar e identical for electrons 
near v = n and composite fermions near v* — n (The 
quantity mw c is the same because it is independent of the 
mass, and £o is t ne same because the magnetic length and 
the disorder correlation lengths are unchanged.), the ra- 
tio of the pinning mode frequency of composite fermions 
and electrons, to®* and w?, respectively, is given by 



cA* 
c*A 



(27) 



Furthermore, cA* and c*A are approximately equal in 
the classical limit, because both c and A are proportional 
to e 2 whereas c* and A* are proportional to e* 2 . There- 
fore, if we take the coupling to the disorder potential 
to be the same, the pinning frequencies of the electron 
WC at v and the CFWC at v* are equal in the classical 
limit, and likely of the same order in the region where 
the system is not too far from the classical region. While 
clearly an approximation, this result is consistent with 
the experimental finding of Zhu et aZ.pand, we believe, 
explains why the pinning frequencies of the electron and 
CFWCs have similar magnitude in a range of parameters. 

We next come to the -B-field dependence of the pin- 
ning frequency. There is an interesting difference in B- 
field dependence of a type-I and type-II WC. For a type-I 
crystal, the B-field dependence of depends on the rel- 
ative magnitudes of I and Tf. The choice — r f gives 
ui® oc 1/B, whereas £o = I produces a B 2 . This dif- 
ference in behavior is easy to identify experimentally and 
the B-field dependence of w° thus provides evidence for 
the relevant length scale for disorder. The B-field depen- 
dence of the type-II WC is more complicated, because 
the B dependence of the CF quasiparticle/quasihole par- 
ticle density p* of the type-II WC must be taken into 
account when examining the _B-field dependence of co®. 
This leads to a simplification for the B dependence of 
the pinning frequency of the CFWC, because the density 
variation with B dominates. Let us consider, for speci- 
ficity, filling factors in the vicinity of v = 1/3,2/5,3/7, 
where v* — n = p<f> /B* and B* — B(l - 2pn). As we 
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move slightly away from the rational fillings, we get: 



= n±v* — 



B* ±SB 



(28) 

where p* = pSB/B* (which is equivalent to the rela- 
tion p* — v* pjv*). For the type-I WC we have (us- 
ing the classical shear modulus) ai® oc p~ 3 ^ 2 , where p 
is independent of B and fixed once the sample is cre- 
ated. For the type-II crystal, when v* < .05, we have 
ujp oc (p*) -3 / 2 oc (SB) ~ 3 / 2 . The dependence of w° on SB 
dominates the _B-ficld dependence of the pinning mode, 
regardless of the choice for the disorder length. This leads 
to the prediction that pinning frequency of the CFWC 
always decreases as one moves away from the rational fill- 
ings, except possibly in the vicinity of the melting tran- 
sition where the shear modulus becomes small. This is 
consistent with the experimental findings^. Similarly, the 
B-field dependence of the pinning modes of the type-II 
crystals observed around integer fillings will also be dom- 
inated by 5B and the frequency of the pinning mode will 
decrease as one moves away from integer fillings, again 
consistent with the experimental observations^'. The 
_B-ficld dependence of the pinning mode thus does not 
shed qualitative light on the nature of disorder for type- 
II WCs. 



D. Melting Temperature 



V* 


Tm 


0.9 


90 mK 


1.06 


25 mK 


1.92 


10 mK 


2.08 


13 mK 


2.94 


2 mK 


3.08 


72 mK 



TABLE I: The peak melting temperatures for CFWC below 
and above v* — 1, 2, and 3, which correspond to electrons 
filling factors in the vicinity of v = 1/3, 2/5, and 3/7. We as- 
sume a density of 1.1 x 10 11 cm -2 and parameters appropriate 
for GaAs. 



Another question of interest is the melting tempera- 
ture of the CFWC. As for two dimensional crystals, it 
is natural to expect that the m elting is described by the 
Kosertlitz-Thouless mechanisir i 40 * 41 !, with melting tem- 
perature given by: 



T M = (2TrV3)- 1 C t 



(29) 



In Fig. [7] we plot the melting temperature of the 
CFWCs near v = 1/3,2/5 and 3/7 in units of e 2 /el m 
50v/5[T] K fa 3.32 x lQ-^y/p/v k for GaAs-AlGaAs 
heterostructures, with magnetic field quoted in Tesla, p 
in cm -2 , and the dielectric function taken to be e = 12.6. 



For a typical density of p = 1.1 x 10 11 cm~ 2 at the fillings 
v = 1/3,2/5,3/7, e 2 /el corresponds to the temperatures 
190K, 174K and 168K, respectively. We also include as 
reference a plot of the melting temperature of the WC of 
point particles with charge e*. For illustration, we give 
in Table I the peak melting temperatures along with the 
filling factors where they occur. 
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point particle WC 
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0.05 0.1 0.15 0.2 
v* 
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. .0002 



CFWC at v* = 3 + v* 


point parti deWC 

1 / 
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v^3/7 



0.05 0.1 0.15 0.2 
— .* 



FIG. 7: Melting temperature Tm as a function of the partial 
filling factor v* — \n — v*\ where v* is the CF filling factor. 
The melting temperature is quoted in units of e 2 /el (with 
fcs = 1). The melting temperature of point particles of charge 
e* is also included for reference. 

Several interesting features are worth noting. The 
melting temperature at 1/3 is in the same ball park as 
the experimental temperature of 50 mK where data in 
Zhu et. al. was collected. However, taken at face value, 
our results would imply a CFWC at the low filling factor 
side of 1/3 but not at filling factors higher than 1/3. The 
reason for the discrepancy is unclear at this stage, and 
may have to do either with some simplifying assumptions 
of the model (such as neglect or disorder, finite width cor- 
rections, or A level mixing), or with slight intrinsic inac- 
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curacies in the estimation of the effective CF interaction. 
We also note that the temperature estimated here is when 
infinitely long range crystalline order is lost; short range 
crystalline order will possibly persist to higher tempera- 
tures. The CFWC melting temperature decreases as we 
go to 2/5 and 3/7, as expected, making their observa- 
tion even more challenging; the reason for anomalously 
high Tm for the CF-hole WC at 3/7 is not understood, 
although it must have to do with the form of the effective 
interaction. An interesting feature is a significant asym- 
metry between the behavior of CFWCs above and below 
v* = integer, as expected from rather general considera- 
tions, given that the interaction between CF particles is 
different from that between CF holes. (Similar asummc- 
try is to be expected for electron WCs above and below 
integer fillings.) The experimental results appear to sug- 
gest a smaller asymmetry than theory; we believe that 
disorder can possibly wash out the asymmetry. 

Also, note that at low fillings, the melting temperature 
for the type-II WC near v = n is (2n + l) 2 larger than 
the corresponding CFWC near v* = n (since c/c* = 
(2n + l) 2 .) The significantly larger melting temperature 
of the type-II WC near v = n explains the relative ease 
with which its pinning mode resonance was detected as 
compared to the resonance near v = 1/3. 



placing a single particle wave packet at each lattice site 
Rj of a hexagonal (or triangular) lattice and antisym- 
metrizing the product. In this section we outline the 
derivation of the effective interaction for the CFWC, 



given in Eq. (12 1. It is not possible to perform this cal- 



culation exactly, and one must resort to approximations. 
We restrict ourselves here to including only the two-body 
direct and exchange terms since these make the greatest 
contribution to the energy.^ 



The expression for the denominator of Eq. ( 11 ) is al- 
ready given to us by |17j : 



where 



<*{K,}({r<})|* { K, } ({r < })> = 



(Al) 



Sij = J rfW^^i^.M^^)^^) 

(A2) 



_IH 2 . 

= e 2 j . 
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Appendix A: Effective model for Maki-Zotos 
Hartree-Fock crystal 



and Rij — \Ri — Rj\. Following appendix A of |17) : 



<* { ^}({n})lE^ ff (r lJ )|vI/ {R3} ({r J }) 



> 



\ £ ^(-Ri, Rj 



(A3) 



As indicated in the text, our calculations are based 
on the Maki-Zotos crystal wave function constructed by 



where 



VzyRii Rj) 



rf 2 r 1 d 2 r 2 ^ tt (r u 0(|V'K I (r-i)nV'H 3 (r-2)| i ~V4 I (^i)^(^)V'H I (r-2)^(r-i)) 



d 2 r 1 d 2 r 2 V eS (r ij )(e-H 



-\[{r 1 -R % ) 2 + (r 2 -R ] ) 2 + (r 1 -R ] ) 2 + {r 2 -R l ) 2 +2i[(x 1 -x 2 )(Y t ~Y ] )-{y 1 -y 2 ){X l -X ] )\]^ 

I 



(A4) 



Transforming to the center of mass and relative coordi- 
nates, the first (direct) term is given by: 



1 



h,ij = ^e—^ I rdrV en {r)e-i r ~I QRijr) (A5) 



,-cir 



1 



and the second (exchange) term is given by: 



= l e ~i R % f rdrV cS (r)e-i r2 J (lRi 3 r) (A6) 



Here, Jo and A) are the Bessel function and the modified 
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Bessel function of the first kind. The nearest neighbor 
term Si 2 is of the order 10~ n for v — .15. Because we 
are not concerned with properties of the system when the 



filling factor is much larger than this, the sum J2i^tj &ij 
is negligible for our purposes. Substituting our results 
for Ii ij and 1% « into Eq. (11) gives: 



<%fl 3 }({r-J)|E^-y eff (%)l^{fl J }({n})> 
<*{*,} ({^})l*{«,}({n})> 



J2i^j h,ij 



h,ij - h,ij 
1 — Sji 



(A7) 



1 - 5« 



±e-i R l frdrV° s (r)e-ir 2 (l (|iV) " Jb(|i2«r)) _ / rdrF eff (r) e -^ (J (|%r) - U\R l3 r)) 



1 - e 



7y R- 



4sinh 



(A8) 



The substitution 1/r for y eff (r) reproduces the Maki- 
Zotos result. For the effective interaction between the 
composite fermions, the integrals can be performed ana- 
lytically, but the expressions are quite cumbersome and 
are most conveniently dealt with in a computer program 
such as Mathematica. 



length), we rewrite this as: 

2 V(R) 



cos(fe-ij)- 



R 



dR x 



d 2 V(R) 27r(n + 771/2) 2irm 



E " am" cos "" i — ~ cos ( B3 ) 



dRl 



Appendix B: Dynamical matrix 

In this section we describe our method for evaluating 
<5 a /3, Eq. (19). We divide $ Qj g into two parts, the k 



independent part and the k dependent part: 



R 



d 2 V{R) 
dR a Rr 



E 

R 



d 2 V(R) 
dR a Rp 



cos(fc--R) (Bl) 



It is necessary to find a way to compute these sums in a 
manner so that the error can be estimated in a reliable 
manner, which is especially important for small k, where 
<fr a p (fe) can be quite small. Since both sums are absolutely 
convergent, we are free to rearrange the terms to produce 
best convergence. 

For specificity, we describe our method for evaluating 
$ XI ; extending this discussion to <fr xy and (& vv is straight- 
forward. We begin with the k dependent term in Eq. 



(Bl ). Our goal is to utilize the periodicity of the cos(fcvR) 



terms to produce a rapidly convergent alternating series. 
Because of inversion symmetry, we have 



cos(fc-i?' 

R 



d 2 V{R) 
dRJ 



where the sum is over the lattice sites R 



^2 cos ( k x R x ) cos ( k y R y ) - 

OrC x 

(B2) 
a (77 + 



^m). Using a = 1 



(with I* as the unit of 



and ^ 



,V3 



k 

2 V 



We 



where 77, m are integers, 4 s = ak 
wish to rearrange the sum on the right hand side of Eq. 
|B3| into an alternating series so that it converges more 
quickly. To this end, we divide the integration region into 
successive light and dark blocks as shown in Fig. [HJ where 
each point represents a lattice vector R. The regions are 

chosen in a manner that cos^ 271 ^"^"^ 2 ^ cos^^ 2 ^ at any 

point in the light (dark) regions gives a positive (nega- 
tive) contribution. Individual positive or negative terms 
in the alternating series are constructed by summing over 
all light or dark blocks along the "diagonal" lines shown 
in Fig. [8] We then apply the Euler transformation^ to 
the resulting alternating series, which produces a rapidly 
converging sum. 

We organize the non fc-dependent sum in a different 
manner. We divide the entire hexagonal lattice into 
successive annuli (labeled by the index 77) as in Fig. |9j 
where, beyond certain distance, the inner radius of each 
successive annulus is twice that of the previous annu- 
lus (e.g. the radii of 5, 10, 20 in Fig. M. The par- 
tial sum over the points in the nth annulus is labeled 

at long distances ensures 



d 2 V{R) 
dRl 



1 

773 



s n . The behavior 
that the sum s n decreases by approximately a factor of 
two for each successive annulus for large 77. Defining 
Sn+i/sn = 1/(2 + e n+ i), with |e n+ i| < |e„|, the sum can 
be expressed exactly as: 



S = 



E = s °+ s i+- • f>+e« + *r (B4) 

R x *=o 



Because we do not know the proper e n a priori, we resort 
to approximating the remaining sum as a power series by 
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n+m/2 



FIG. 8: A portion of the upper right hand quadrant of points 
used in the summation of Eq. ( B3 1 . The points on the graph 
mark the locations (n + m/2, m), where n, m are positive in- 
tegers that occur in the expression Eq. (JoJ . The quantities 
t x and t y are defined in Appendix [B] 



setting e 



0. The sum is then given by: 



S' 



d 2 V(R) 

R 



dm 



S0+SlH \S; 



E 2 " 

i=0 



So+Sl 



2s n 



(B5) 

To estimate the error in this approximation, we define 
S ± = s + ai + • • • + s„ E^o( 2 T knl)"' so that S~ < 
S < S + . With some algebra, and assuming e„ <?C 1, this 
inequality becomes: 




FIG. 9: Concentric annuli illustrate how the sum in Eq. (B4| 
is organized. The points on the graph have the same locations 
as in Fig. [8] 



s„(2-|e„|) < 5-s -si s„_i < s„(2+|e„|), (B6) 



which establishes the error to be bounded by \S — S'\ < 
Sn\£n\ — |2s„ — s„_i| In our numerical calculations, e„ 
becomes rapidly smaller as we increase n, allowing us to 
calculate the sum to a very high degree of accuracy. We 
implement the sum using the epsilon algorithm given in 
Ref. 02]. 
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